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TECHNIQUES FOR IDENTIFYING, CONFIRMING, MAPPING, AND 
5 CATEGORIZING POLYMERS 

This application claims the benefit of U.S. Provisional Application No. 60/057,957, 
filed September 5, 1997, which is hereby incorporated by reference. 

1 0 GOVERNMENT RIGHTS NOTICE 

Portions of the material in this specification arose under the cooperative agreement 
70NANABSH1031 between Affymetrix, Inc. and the Department of Commerce through the 
Advanced Technology Program. 



1 5 BACKGROUND OF THE INVENTION 

The present invention relates to detecting differences in polymers. More specifically, 
the present invention relates to techniques for identifying, confirming, mapping, and 
genotyping sample nucleic acid sequences. 

Devices and computer systems for forming and using arrays of materials on a chip or 
20 substrate are known. For example, PCT applications W092/1 0588 and 95/1 1 995, both 
incorporated herein by reference for all purposes, describe techniques for sequencing or 
sequence checking nucleic acids and other materials. Arrays for performing these operations 
may be formed according to the methods of, for example, the pioneering techniques disclosed 
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in U.S. Patent Nos. 5,445,934, 5,384,261 and 5,571,639, each incorporated herein by 
reference for all purposes. 

According to one aspect of the techniques described therein, an array of nucleic acid 
probes is fabricated at known locations on a chip. A labeled nucleic acid is then brought into 
5 contact with the chip and a scanner generates an image file indicating the locations where the 
labeled nucleic acids are bound to the chip. Based upon the image file and identities of the 
probes at specific locations, it becomes possible to extract information such as the nucleotide 
or monomer sequence of DNA or RNA. Such systems have been used to form, for example, 
arrays of DNA that may be used to study and detect mutations relevant to genetic diseases, 
1 0 cancers, infectious diseases, HIV, and other genetic characteristics. 
! The VLSIPS™ technology provides methods of making very large arrays of 

! oligonucleotide probes on very small chips. See U.S. Patent No. 5,143,854 and PCT patent 

| publication Nos. WO 90/15070 and 92/10092, each of which is incorporated by reference for 

all purposes. The oligonucleotide probes on the DNA probe array are used to detect 

| 1 5 complementary nucleic acid sequences in a sample nucleic acid of interest (the 'target" 

i 
I 

nucleic acid). 

For sequence checking applications, the chip may be tiled for a specific target nucleic 
acid sequence. As an example, the chip may contain probes that are perfectly complementary 
to the target sequence and probes that differ from the target sequence by a single base 
20 mismatch. For de novo sequencing applications, the chip may include all the possible probes 
of a specific length. The probes are tiled on a chip in rows and columns of cells, where each 
cell includes multiple copies of a particular probe. Additionally, "blank" cells may be present 
on the chip which do not include any probes. As the blank cells contain no probes, labeled 
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targets should not bind specifically to the chip in this area. Thus, a blank cell provides a 
measure of the background intensity. 

While the Human Genome Project is attempting to produce the first complete 
reference sequence of the human chromosomes, attention is already focusing on the sequence 
5 variations among individuals. The genetic diversity is of interest because it may explain the 
basis of heritable variation in disease susceptibility, as well as harboring a record of human 
genetic migrations. 

The most common type of human genetic variation is the single-nucleotide 
polymorphism (SNP), which is a position where two alternative bases occur at appreciable 
1 0 frequency (e.g. , greater than 1 %) in the human population. There are many uses for SNPs 
including serving as genetic markers for identifying disease genes by linkage studies in 
families, linkage disequilibrium in isolated populations, association analysis of patients and 
controls, and loss-of-heterozygosity studies in tumors to name a few. It is believed that large 
collections of mapped SNPs would provide a powerful tool for human genetic studies. 
15 Although individual SNPs can be less informative than conventional genetic markers, SNPs 
can be more abundant and have a greater potential for automation. 

Accordingly, there is a need for innovative techniques for identifying, confirming, 
mapping, and categorizing polymers, such as nucleic acids. 
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SUMMARY OF THE INVENTION 
The present invention provides innovative systems and methods for identifying, 
confirming, mapping, and categorizing polymers, such as nucleic acids. Samples of polymers 
5 can be analyzed to determine estimations of the fractions of distinct polymers in the samples. 
The estimations can then be clustered in order to form a fraction pattern. The fraction pattern 
can be utilized to interpret data from subsequent samples in order to more accurately identify, 
confirm, map, and/or categorize the polymers in these samples. Several embodiments of the 
invention are described below. 

10 In one embodiment, the invention provides a method of estimating a fraction of first 

and second polymers in a sample of polymers. A first hybridization value indicative of 
hybridization affinity of the sample of polymers to polymer probes that are complementary to 
the first polymer is input. A second hybridization value indicative of hybridization affinity of 
the sample of polymers to polymer probes that are complementary to the second polymer is 

15 input. The estimation of the fraction the first and second polymers in the sample of polymers 
is calculated by dividing the first hybridization value by a sum of the first and second 
hybridization values. 

In another embodiment, the invention provides a method of estimating a fraction of 
alleles in a sample of nucleic acids. A first hybridization value indicative of hybridization 
20 affinity of the sample of nucleic acids to nucleic acid probes that are complementary to a 
portion of a first allele is input. A second hybridization value indicative of hybridization 
affinity of the sample of nucleic acids to nucleic acid probes that are complementary to a 
portion of a second allele is input. The estimation of the fraction of the first and second 
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alleles in the sample of nucleic acids is calculated by dividing the first hybridization value by 
a sum of the first and second hybridization values. 

In another embodiment, the invention provides a method of genotyping a sample of 
nucleic acids. Multiple estimations of the fraction of first and second alleles in multiple 
samples of nucleic acids are generated. The multiple estimations are clustered to form a 
fraction pattern. An estimation of the fraction of the first and second alleles in a sample of 
nucleic acids is calculated and the sample of nucleic acids is genotyped by comparing the 
estimation of the fraction of the first and second alleles in the sample of nucleic acids to the 
fraction pattern. As an example, the clusters in the fraction pattern can represent homozygous 
for a first allele, homozygous for a second allele or heterozygous for both alleles. 

Other features and advantages of the invention will become readily apparent upon 
review of the following detailed description in association with the accompanying drawings. 
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BRIEF DESCRIPTION OF THE DRAWINGS 
Fig. 1 illustrates an example of a computer system that may be utilized to execute the 
software of an embodiment of the invention. 
5 Fig. 2 illustrates a system block diagram of the computer system of Fig. 1 . 

Fig. 3 illustrates an overall system for forming and analyzing arrays of biological 
materials such as DNA or RNA. 

Fig. 4 illustrates conceptually the binding of probes on chips. 
Fig. 5 illustrates a portion of one embodiment of a genotyping chip. 
10 Fig. 6 shows a flowchart of a process of estimating fractions of polymers. 

Fig. 7 shows a flowchart of a process of analyzing a sample of nucleic acids. 
Fig. 8 is a table of experimental data and results from one embodiment of the 
invention. 

Fig. 9A-9C illustrate fraction patterns from the table of Fig. 8. 
15 Fig. 1 0 shows a flowchart of a process of making a genotype call form a single strand. 

Fig. 1 1 shows a flowchart of a process of combining strand genotype calls into a single 
genotype call. 
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DETAILED DESCRIPTION OF PREFERRED EMBODIMENTS 

In the description that follows, the present invention will be described in reference to 

5 preferred embodiments that utilize VLSIPS™ technology for making very large arrays of 
oligonucleotide probes on chips. However, the invention is not limited to nucleic acids or to 
this technology and may be advantageously applied to other polymers and manufacturing 
processes. Therefore, the description of the embodiments that follows for purposes of 
illustration and not limitation. 

10 Fig. 1 illustrates an example of a computer system that may be used to execute the 

software of an embodiment of the invention. Fig. 1 shows a computer system 1 that includes 
a display 3, screen 5, cabinet 7, keyboard 9, and mouse 1 1 . Mouse 1 1 may have one or more 
buttons for interacting with a graphical user interface. Cabinet 7 houses a CD-ROM drive 13, 
system memory and a hard drive (see Fig. 2) which may be utilized to store and retrieve 

1 5 software programs incorporating computer code that implements the invention, data for use 
with the invention, and the like. Although a CD-ROM 1 5 is shown as an exemplary computer 
readable storage medium, other computer readable storage media including floppy disk, tape, 
flash memory, system memory, and hard drive may be utilized. Additionally, a data signal 
embodied in a carrier wave (e.g. , in a network including the Internet) may be the computer 

20 readable storage medium. 

Fig. 2 shows a system block diagram of computer system 1 used to execute the 
software of an embodiment of the invention. As in Fig. 1, computer system 1 includes 
monitor 3 and keyboard 9, and mouse 11. Computer system 1 further includes subsystems 
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such as a central processor 51, system memory 53, fixed storage 55 (e.g., hard drive), 
removable storage 57 (e.g., CD-ROM drive), display adapter 59, sound card 61, speakers 63, 
and network interface 65. Other computer systems suitable for use with the invention may 
include additional or fewer subsystems. For example, another computer system could include 
5 more than one processor 5 1 (/. e. , a multi-processor system) or a cache memory. 

The system bus architecture of computer system 1 is represented by arrows 67. 
However, these arrows are illustrative of any interconnection scheme serving to link the 
subsystems. For example, a local bus could be utilized to connect the central processor to the 
system memory and display adapter. Computer system 1 shown in Fig. 2 is but an example of 
10 a computer system suitable for use with the invention. Other computer architectures having 
different configurations of subsystems may also be utilized. 

For purposes of illustration, the present invention is described as being part of a 
computer system that designs a chip mask, synthesizes the probes on the chip, labels the 
nucleic acids, and scans the hybridized nucleic acid probes. Such a system is fully described 
15 in U.S. Patent No. 5,571 ,639 that has been incorporated by reference for all purposes. 

However, the present invention may be used separately from the overall system for analyzing 
data generated by such systems. 

Fig. 3 illustrates a computerized system for forming and analyzing arrays of biological 
materials such as RNA or DNA. A computer 100 is used to design arrays of biological 
20 polymers such as RNA and DNA. The computer 100 may be, for example, an appropriately 
programmed Sun Workstation or personal computer or workstation, such as an IBM PC 
equivalent, including appropriate memory and a CPU as shown in Figs. 1 and 2. The 
computer system 100 obtains inputs from a user regarding characteristics of a gene of interest, 
and other inputs regarding the desired features of the array. Optionally, the computer system 
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may obtain information regarding a specific genetic sequence of interest from an external or 
internal database 102 such as GenBank. The output of the computer system 100 is a set of 
chip design computer files 104 in the form of, for example, a switch matrix, as described in 
PCT application WO 92/10092, and other associated computer files. 
5 The chip design files are provided to a system 1 06 that designs the lithographic masks 

used in the fabrication of arrays of molecules such as DNA. The system or process 106 may 
include the hardware necessary to manufacture masks 1 1 0 and also the necessary computer 
hardware and software 108 necessary to lay the mask patterns out on the mask in an efficient 
manner. As with the other features in Fig. 3, such equipment may or may not be located at the 

10 same physical site but is shown together for ease of illustration in Fig. 3. The system 106 
generates masks 1 10 or other synthesis patterns such as chrome-on-glass masks for use in the 
fabrication of polymer arrays. 

The masks 1 10, as well as selected information relating to the design of the chips from 
system 100, are used in a synthesis system 112. Synthesis system 1 12 includes the necessary 

1 5 hardware and software used to fabricate arrays of polymers on a substrate or chip 1 14. For 
example, synthesizer 1 12 includes a light source 1 1 6 and a chemical flow cell 1 1 8 on which 
the substrate or chip 1 14 is placed. Mask 1 10 is placed between the light source and the 
substrate/chip, and the two are translated relative to each other at appropriate times for 
deprotection of selected regions of the chip. Selected chemical regents are directed through 

20 flow cell 1 1 8 for coupling to deprotected regions, as well as for washing and other operations. 
All operations are preferably directed by an appropriately programmed computer 119, which 
may or may not be the same computer as the computer(s) used in mask design and mask 
making. 
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The substrates fabricated by synthesis system 1 12 are optionally diced into smaller 
chips and exposed to marked targets. The targets may or may not be complementary to one or 
more of the molecules on the substrate. The targets are marked with a label such as a 
fluorescein label (indicated by an asterisk in Fig. 3) and placed in scanning system 120. 
5 Although preferred embodiments utilize fluorescent markers, other markers may be utilized 
that provide differences in radioactive intensity, light scattering, refractive index, 
conductivity, electroluminescence, or other large molecule detection data. Therefore, the 
present invention is not limited to analyzing fluorescence measurements of hybridization but 
may be readily utilized to analyze other measurements of hybridization. 

1 0 Scanning system 120 again operates under the direction of an appropriately 

programmed digital computer 122, which also may or may not be the same computer as the 
computers used in synthesis, mask making, and mask design. The scanner 120 includes a 
detection device 124 such as a confocal microscope or CCD (charge-coupled device) that is 
used to detect the location where labeled target (*) has bound to the substrate. The output of 

1 5 scanner 120 is an image file(s) 124 indicating, in the case of fluorescein labeled target, the 
fluorescence intensity (photon counts or other related measurements, such as voltage) as a 
function of position on the substrate. Since higher photon counts will be observed where the 
labeled target has bound more strongly to the array of polymers (e.g., DNA probes on the 
substrate), and since the monomer sequence of the polymers on the substrate is known as a 

20 function of position, it becomes possible to determine the sequence(s) of polymer(s) on the 
substrate that are complementary to the target. 

The image file 124 is provided as input to an analysis system 126 that incorporates the 
synthesis integrity evaluation techniques of the present invention. Again, the analysis system 
may be any one of a wide variety of computer system(s), but in a preferred embodiment the 

10 



WO 99/14369 



PCT/US98/18873 



analysis system is based on a WINDOWS NT workstation or equivalent. The analysis system 
may analyze the image file(s) to generate appropriate output 128, such as the identity of 
specific mutations in a target such as DNA or RNA. 

Fig. 4 illustrates the binding of a particular target DNA to an array of DNA probes 
5 1 1 4. As shown in this simple example, the following probes are formed in the array: 
3 ' -AGAACGT 
AGACCGT 
AGAGCGT 
AGATCGT 
10 o 

o 
o 

As shown, when the fluorescein-labeled (or otherwise marked) target 5 f -TCTTGCA is 
exposed to the array, it is complementary only to the probe 3 1 -AGAACGT, and fluorescein 
15 will be primarily found on the surface of the chip where 3 ' -AGAACGT is located. The chip 
contains cells that include multiple copies of a particular probe and the cells may be square 
regions on the chip. 

In the description that follows, the invention will be described in reference to SNPs. 
Therefore, at these SNP positions, the sample nucleic acids can generally include one or both 
20 of two alleles, which will be referred to generically as the A-allele and the B-allele. In a pool 
of individuals, it would be expected that if there is an SNP at a location, some individuals 
would be homozygous A A, some would be homozygous BB and some would be heterozygous 
AB. Although the invention will be described in reference to SNPs, the invention is not so 
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limited and can be advantageously applied to other polymorphism conditions including 
insertions, deletions, multi-position polymorphisms, multi-base polymorphisms, and multi- 
allelic positions. 

Fig. 5 shows a portion of a tiled chip that can be utilized to identify, confirm, map, 
5 and/or genotype sample nucleic acids. The chip typically includes thousands of cells where 
nucleic acid probes are coupled to the chip. The chip can have different probes tiled on the 
chip for a specific applications (e.g., detecting mutations that are indicative of a change in 
resistance to a drug therapy in treating HIV) or the chips can be generic (e.g., have all possible 
8-mer probes). Furthermore, the cells can be laid out on the chip in any number of ways. For 

10 example, the cells can be tiled in locations to improve scan data by reducing inter-cell 

interference, to reduce the cost of manufacturing by decreasing the number of masks utilized, 
or to make it easier for users to interpret the data. 

A preferred layout of a portion of genotyping chip will be described in reference to 
Fig. 5. The layout is generally organized so it is easier for users to interpret the data and 

1 5 understand the structure of the probes on the chip as will be described herein. However, it 
should be understood that the invention is not limited to any specific chip layout or structure 
of the probes. The invention can be advantageously applied to other chip layouts as well. 

Assume that the A-allele includes the sequence S'-CCTCGGCAT and the B-allele 
includes the sequence S'-CCTCAGCAT, where the underlined nucleotide represents the SNP 

20 that distinguishes the alleles. For discussion purposes, the position of nucleotides in the 

alleles will be specified by the variable d, where d = 0 for the SNP and d decreases to the left 
of the SNP and increases to the right of the SNP. Fig. 5 shows a block 201 that includes five 
positions around the SNP. The number of positions included in block 201 has been chosen 
for illustration purposes and can be varied depending on the application. 

12 
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Along the top of Fig. 5 is shown the sequence of the alleles. At d = 0, the nucleotide is 
shown as G/A to represent the SNP that distinguishes the A-allele and the B-allele. As 
shown, block 201 includes four rows of cells tiled for the A-allele and four rows of cells tiled 
for the B-allele. The rows of cells are designed by nucleotides T, G, A, and A, which 
5 represent the nucleotide at the interrogation position of the nucleic acid probes in that row. 
Block 201 can be further divided into miniblocks 203, which represent a set of cells (and their 
corresponding probes) that have a perfect match probe for the alleles of interest and zero or 
more mismatch probes. Although the miniblocks are shown as columns on the chip in one 
embodiment, they are not required to be in any particular spatial orientation on the chip. 
1 0 Recalling that the sequence of the alleles is 5'-CCTC?GCAT where the question mark 

indicates that the SNP that can be either a G or A. Assuming that the probes on the chip are 
5-mers with an interrogation position at the center nucleotide, the probes on the chip would be 
as follows: 





d 


-2 


-1 


0 


1 


2 


15 


A-allele 


GGTGC 


GAJCC 


AGTCG 


GCTGT 


CCTTA 






GGGGC 


GAGCC 


AGGCG 


GCGGT 


CCTTA 






GGCGC 


GACCC 


AGCCG 


GCCGT 


CCTTA 






GGAGC 


GAACC 


AGACG 


GCAGT 


CCTTA 


20 


B-allele 


GGTGT 


GATTC 


AGTCG 


GTTGT 


TCTTA 






GGGGT 


GAGTC 


AGGCG 


GTGGT 


TCTTA 






GGCGT 


GACTC 


AGCCG 


GTCGT 


TCTTA 






GGAGT 


GAATC 


AGACG 


GTAGT 


TCTTA 



25 The interrogation position nucleotides are underlined in the probes. Additionally, the SNP 
nucleotides are bolded in each probe : C in the case of the A-allele and T in the case of the B- 
allele. As shown, the miniblock at d = 0 has two identical halves, accordingly it would be 
possible to only tile one half of the miniblock at that location. 
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The miniblocks shown include a perfect match (i.e., a probe perfectly complementary 
to the target allele) to the target alleles and multiple probes with a single mismatch. The 
center miniblocks are unique in that the interrogation position is the same as the SNP. 
Accordingly, the mismatch probes of the center miniblock have a single base mismatch from 
5 both of the alleles. This is in contrast to the noncenter miniblocks where the mismatch probes 
have a single base mismatch to one allele and two base mismatches from the other allele. 

Returning to Fig. 5, the cells are labeled with the following: 

PMa - perfect match to A-allele 

PMb - perfect match to B-allele 
10 Ma - single base mismatch to A-allele 

Mb - single base mismatch to B-allele 

M* - single base mismatch to both the A-allele and B-allele 
Fig. 5 illustrates that the noncenter miniblocks include two perfect match probes and four 
mismatch probes, two for each allele. The center miniblocks include perfect match probes for 
1 5 each allele and two mismatch probes for both alleles. In the discussion that follows, these 
labels will be utilized to describe embodiments of the invention. 

Fig. 6 is a flowchart of a estimating fractions of polymers in a sample having multiple 
polymers. Although in preferred embodiments, the polymers are nucleic acids, other 
polymers including peptides, oligosaccharides and the like can be utilized. Therefore, the 
20 invention is not limited to the specific embodiments described herein. 

At a step 301, a first hybridization value indicative of hybridization affinity of a 
sample of polymers to a first polymer is input. The sample of polymers can have multiple 
different sequences of polymers including those that are complementary to the first polymer. 
A second hybridization value indicative of hybridization affinity of a sample of polymers to a 
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first polymer is input at a step 303. The first and second hybridization values are expected to 
vary according to the fraction of the sample of polymers that are complementary to the first 
and second polymer. 

At a step 305, a third hybridization value is subtracted from the first and second 
5 hybridization values. The third hybridization value can be utilized to normalize the first and 
second hybridization values. As will be described below, in preferred embodiments, the third 
hybridization value is derived from mismatched nucleic acid probes, which can make 
subtracting a background value unnecessary. 

An estimation of the fraction of the first and second polymers in the sample of 

10 polymers is calculated by dividing the first normalized hybridization value by the sum of the 
first and second normalized hybridization values. In this manner, the estimation will 
generally vary between 0 and 1, where 1 indicates a very high concentration of polymers 
complementary to the first polymer and 0 indicates a very high concentration of polymers 
complementary to the second polymer. 

15 In embodiments where more than one polymer or allele is being analyzed, the 

denominator can include three or more hybridization values. Accordingly, the invention is not 
limited to analyzing only two polymers or alleles. 

Returning back to the labels of the nucleic acid probes of a preferred embodiment, it 
would be expected that if a sample of nucleic acid probes was taken from an individual that 

20 was homozygous for the A-allele, then PMa > PMb = M* at the center miniblock. Similarly, 
it would be expected that if a sample of nucleic acid probes was taken from an individual that 
was homozygous for the B-allele, then PMb > PMa = M* at the center miniblock. If the 
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individual was heterozygous, it would be expected that PMa = PMb > M* at the center 
miniblock. 

The estimation of the fraction shown in Fig. 6 can be as follows for each of the 
miniblocks: 

5 p-hat = (PMa - MM) / (PMa - MM) + (PMb - MM) 

where p-hat represents the estimation of the fraction and MM is derived from the mismatch 
probes M*, Ma or Mb. As an example, MM can be calculated as follows: 
MM = min(M*, M*) for the center miniblocks 

MM = max(median(Ma, Ma, Ma), median(Mb, Mb, Mb)) for the noncenter 
10 miniblocks 

The minimum is utilized for center miniblocks in one embodiment because there are only two 

distinct M* mismatch probes and it can be desirable to avoid including a very large value that 

could result from cross-hybridization events. 

The above equations for MM are illustrative of those that can be utilized. Others can 
15 be utilized and may be preferably varied depending on the specific application. As an 

example, the miniblocks can be tiled with only one mismatch probe. In such a case, MM can 

be equal to the hybridization affinity measured from probes in the single mismatch cell. 

Although mismatch probes are not required, they have been found to make the calculations 

more robust. 

20 It will be apparent from looking at the above equation for p-hat that if PMa or PMb is 

less than MM, p-hat can become negative and possibly result in a division by zero. 
Accordingly, in some embodiments, if a calculation of (PMa - MM) or (PMb - MM) results in 
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a negative number, the result is set equal to zero. Accordingly, if Da = max(PMa - MM, 0) 
and Db - max(PMb - MM, 0), the above p-hat formula can be written as follows: 
p-hat = Da/(Da + Db) 

It would still be possible to have a denominator of zero, but as will be described in more detail 
5 later, a quality filter can reject such data as bad so the illegal operation does not occur. 
Another alternative would be to set a negative (or zero) result to a very small number like 
0.00001, as an example. 

In another embodiment, the fraction estimation can be calculated from multiple 
miniblocks. For example, the hybridization affinity of the perfect match probes PMa from the 

10 miniblocks can be summed and the hybridization affinity of the perfect match probes PMb 
from the miniblocks can be summed. The sums are then utilized in one of the above formulas 
as PMa and PMb, respectively, to produce the fraction estimation p-hat. One advantage of 
this embodiment is that the miniblocks are effectively weighted according to the hybridization 
affinities of the perfect match probes, which can be desirable in some applications. 

1 5 Fig. 7 shows a flowchart of a process of analyzing a sample of multiple nucleic acids. 

At a step 351, quality tests are performed on the hybridization affinity data. As in all the 
flowcharts herein, there are many steps that can be deleted, added, combined, and modified 
while still being within the spirit of the invention. Quality tests are not required in all 
embodiments but they can be beneficial in removing hybridization affinity data that are 

20 erroneous (e.g., there is a flaw in the chip or an error occurred while scanning a region of the 
hybridized chip). One quality test that can be utilized is that (PMa - MM) and (PMb - MM) 
must be above some predetermined threshold (e.g., 15 photon counts). Another quality test 
can be the highest hybridization affinity must be from a perfect match probe PMa or PMb. 
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Another quality test can be if (PMa / max(Ma) < 1 and (PMb / max(Mb) < 1 , then the 
miniblock fails. Other quality tests can be devised and utilized. Additionally, the quality tests 
can be combined so that more than one must be passed for further processing. Quality tests 
can also be performed on blocks as a whole if desired. 
5 Multiple estimations of the fraction of alleles in samples are generated at a step 353. 

Each estimation of the fraction of nucleic acids in a sample can be calculated according to the 
process shown in Fig. 6 or described above. The samples can be taken from a reasonably 
large and diverse group of individuals. The multiple estimations of the fraction of alleles in 
the samples are clustered to form a fraction pattern at a step 355. The clustering can be 

1 0 performing in a conventional maimer including the valley seeking clustering techniques 
described by W.L.G. Koontz et al. IEEE Trans., Comp. C-21, 171 (1972), which is hereby 
incorporated by reference. Other clustering techniques can be utilized with the invention. As 
the following will describe, the fraction pattern can be subsequently utilized for identifying, 
confirming, mapping, and/or genotyping samples of nucleic acids. 

15 At a step 357, a sample is analyzed and an estimation of the fraction of alleles in the 

sample is calculated. The estimation can be calculated according to the process shown in Fig. 
6 and described above. The estimation of the fraction of alleles is compared to the fraction 
pattern at a step 359. The comparison can result in identifying an SNP, confirming its 
existence, mapping the target, and/or genotyping the sample (or individual from which it was 

20 obtained). 

It may be beneficial to now describe data from an experiment performed with an 
embodiment of the invention. A genotyping chip having a structure similar to that described 
in reference to Fig. 5 was hybridized to samples from 41 individuals. The chip included 
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approximately 600 SNPs. Fig. 8 shows a tabular summary of the results of the experiment for 
the first 45 SNPs. 

A table 375 includes multiple columns of which column 376 is the number of the SNP 
(here 1-45) that is depicted in the designated row of the table. Column 377 is the name of the 
5 SNP (or marker). Although to reduce confusion, the nucleic acids have been discussed as 
though there was only one direction, preferred embodiments utilize sample nucleic acids and 
probes of both the sense and antisense strands. Column 378 specifies from which strand the 
data represents with "A" representing sense and "Z" representing antisense. 

Column 379 depicts the number of samples that were not utilized for the particular 
10 SNP (e.g., the data did not pass a quality filter). Column 380 indicates the number of 

individuals that are homozygous AA, column 381 indicates the number of individuals that are 
heterozygous AB, and column 382 indicates the number of individuals that are homozygous 
BB. Column 383 indicates the center of the clusters in the fraction pattern that corresponds to 
the homozygous AA individuals, column 384 indicates the center of the clusters in the 
1 5 fraction pattern that corresponds to the heterozygous AB individuals, and column 385 

indicates the center of the clusters in the fraction pattern that corresponds to the homozygous 
BB individuals. The estimations of the fractions were multiplied by 1 00 to give values for the 
fraction pattern that range from approximately 0 to 100 as shown for convenience. 

It may be expected that the center of the clusters for homozygous AA individuals 
20 would be at 100, the center of the clusters for the heterozygous AB individuals would be at 
50, and the center of the clusters for homozygous BB individuals would be at 0. A perusal of 
table 375 shows that this is rarely the case with real data. A row 391 shows very good data in 
that the center of the clusters for columns 383-385 is at 99, 48 and 1, which is very close to 
what one might expect. However, a row 393 includes data where the center of the clusters for 
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columns 383-385 is at 97, 71 and 10, which is very different from what one would probably 
expect. 

Data like that shown in row 393 illustrates that it may not be entirely accurate to 
genotype according to a strict 100, 50 and 0 scale. The fraction pattern may be unique for 
5 each SNP. For example, the fraction pattern may be the result of the specific nucleic acid 
sequences near the SNP, the interaction of the nucleic acids during hybridization and the like. 
In order to genotype a sample (or individual) at the SNP corresponding to row 393, 
embodiments of the invention utilize the fraction pattern of 97, 71 and 10. Accordingly, if the 
estimation of the fraction of alleles in the sample is within a predefined range of 71, the 
10 sample is genotyped as heterozygous. Thus, an estimation of 75 could result in a genotype 
call of heterozygous, whereas using a 100, 50 and 0 scale may likely result in noncall or an 
erroneous call. 

The existence and relative positions of three clusters can be utilized to identify and/or 
confirm SNPs. SNPs can also be identified and confirmed according to the innovative 

1 5 techniques describe in U.S. Patent Application No. , filed August 14, 1 998, entitled 

"Polymorphism Detection Utilizing Clustering Analysis," Attorney Docket No. AFFYP010, 
which is hereby incorporated by reference. Other techniques can be utilized to identify and 
confirm SNPs including conventional gel-based techniques. 

A row 395 shows only two clusters, as there is no cluster for homozygous BB. This 

20 could be the result of not having a sample from an individual that is homozygous BB. 

However, unusual clustering patterns can also demonstrate that there is not an SNP at that 
position or that a more complex mutation is occurring. For example, a single cluster can 
indicate that there is not an SNP at that position, whereas a clustering of more than three 
cluster can be indicative of more differences in the nucleic acids than a single polymorphism. 
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Figs. 9A-9C illustrate the fraction patterns of rows 391, 393 and 395, respectively. As 
shown, Fig. 9A looks to be almost ideal (z.e., expected) in that the clusters are well spaced 
apart. However, Fig. 9B shows a fraction pattern that has two clusters fairly close together. 
And lastly, Fig. 9C shows a fraction pattern that only has two clusters. The sizes of the 
5 clusters do not depict the relative number of samples within the cluster but are shown to 
illustrate to the reader how different the center of the clusters can be and how the fraction 
patterns can vary depending on the SNP. 

The relative change in fraction patterns for an individual can be monitored over time. 
For example, a tumor causing SNP can be monitored to determine if the relative 

1 0 concentrations of the allele believed responsible for the tumor is decreasing with therapy. 

Fig. 1 0 shows a flowchart of a process of making a single strand genotype call. At a 
step 401 , a check can be performed to determine if two or more miniblocks had bad data. For 
example, two or more miniblocks for the SNP did not pass a quality test of step 35 1 of Fig. 7. 
If there was two or more miniblocks with bad data at a step 403, the genotype call is "fail" at a 

15 step 405. Otherwise, a check can be performed to determine if the most common genotype 
call from the miniblocks passes a test or tests at a step 407. Although the check in step 401 
has been described in reference to two or more miniblocks, this number can be varied in 
depending on the application. However, it may be beneficial if the number of miniblocks for 
the test is varied proportionally with the number of miniblocks (i.e. 9 if the number of 

20 miniblocks increases, so does the number involved in the check). 

A test be n > (N - s) / 2, where N is equal to the number of miniblocks, n is the number ' 

i 
i 

of miniblocks making the most common genotype call and s is the number of miniblocks that 
had bad data or did not pass a quality test. Another test can be m < 2, where m is equal to the 
number of miniblocks making the second most common genotype call. Additionally, a test 
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can be n - m > (N - s) / 2. If the miniblocks fail to pass one or a combination of the tests (e.g., 
all the equations are true) at a step 409, the genotype call is "no call" at a step 411. Other tests 
and combinations of tests can also be devised and utilized. If the miniblocks pass the test or 
tests, the genotype call is the most common genotype call made by the miniblocks at a step 
5 413. 

Fig. 1 1 shows a flowchart of a process of combining strand genotype calls. The 
multiple strand genotype calls can be made by a process as shown in Fig. 10. For simplicity, 
the flowchart of Fig. 1 1 will be described in reference to two strands. However, the process 
can be advantageously applied to combining genotype calls from more than two strands and to 

1 0 results from different experiments. 

At a step 451, a check can be performed to determine if both strands produced a "no 
call." If both strands did produce a "no call" at a step 453, the combined call is a "no call" at 
a step 455. Otherwise, a check is performed to determine if one strand produced a "no call" 
and the other produced a genotype call (e.g., homozygous AA). If the strands had a mixed 

1 5 "no call" result, the genotype call is made at a step 46 1 . 

A check is performed at a step 463 to determine if both strands make a genotype call 
and agree on the call. If both strands do not produce the same genotype call, the call is a "no 
call" at a step 467. Otherwise, the genotype call for both strands is made as the agreed 
genotype call at a step 469. 

20 Embodiments of the invention have been utilized to simultaneously genotype 500 

SNPs utilizing a single genotyping chip. See "Large-Scale Identification, Mapping, and 
Genotyping of Single-Nucleotide Polymorphisms in the Human Genome," David G. Wang et 
al., Science, Vol. 280, pp. 1077-1082, May 15, 1998, which is hereby incorporated by 
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reference. The results demonstrated the feasibility of large-scale identification and analysis of 
SNPs was possible utilizing embodiments of the invention. 

While the above is a complete description of preferred embodiments of the invention, 
various alternatives, modifications, and equivalents may be used. It should be evident that the 
invention is equally applicable by making appropriate modifications to the embodiments 
described above. For example, the invention has been described in reference to nucleic acid 
probes that are synthesized on a chip. However, the invention may be advantageously applied 
to other monomers (e.g., amino acids and saccharides) and other hybridization techniques 
including those where the probes are not attached to a substrate. Therefore, the above 
description should not be taken as limiting the scope of the invention that is defined by the 
metes and bounds of the appended claims along with their full scope of equivalents. 
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CLAIMS 

What is claimed is: 

1 . A method of estimating a fraction of first and second polymers in a sample of 
5 polymers, comprising: 

inputting a first hybridization value indicative of hybridization affinity of the sample 
of polymers to polymers probes that are complementary to the first polymer; 

inputting a second hybridization value indicative of hybridization affinity of the 
sample of polymers to polymers probes that are complementary to the second polymer; and 
1 0 calculating an estimation of the fraction the first and second polymers in the sample of 

polymers by dividing the first hybridization value by a sum of the first and second 
hybridization values. 

2. The method of claim 1 , further comprising subtracting from the first and 

1 5 second hybridization values a third hybridization value indicative of hybridization affinity of 
the sample of polymers to polymers probes that include at least one monomer that differs from 
both the first and second polymers. 

3. The method of claim 2, further comprising selecting the third hybridization 
20 value from a plurality of hybridization values, each of which is indicative of hybridization 

affinity of the sample of polymers to polymers probes that include at least one monomer that 
differs from both the first and second polymers. 

4. The method of claim 3, wherein selecting the third hybridization value includes 
25 selecting a hybridization value indicative of lowest hybridization affinity to the sample of 

polymers. 

5. The method of claim 2, wherein the subtracted first and second hybridization 
values are set equal to zero if they are less then zero. 

30 
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6. The method of claim 1 , further comprising: 

generating a plurality of estimations of the fraction the first and second polymers in a 
plurality of samples of polymers; and 

clustering the plurality of estimations to form a fraction pattern. 

7. A computer program product that estimates a fraction of first and second 
polymers in a sample of polymers, comprising: 

computer code that receives a first hybridization value indicative of hybridization 
affinity of the sample of polymers to polymers probes that are complementary to the first 
polymer; 

computer code that receives a second hybridization value indicative of hybridization 
affinity of the sample of polymers to polymers probes that are complementary to the second 
polymer; 

computer code that calculates an estimation of the fraction of the first and second 
polymers in the sample of polymers by dividing the first hybridization value by a sum of the 
first and second hybridization values; and 

a computer readable medium that stores the computer codes. 

8. The computer program product of claim 7, wherein the computer readable 
medium is selected from the group consisting of floppy disk, tape, flash memory, system 
memory, hard drive, and a data signal embodied in a carrier wave. 

9. A method of estimating a fraction of alleles in a sample of nucleic acids, 
comprising: 

inputting a first hybridization value indicative of hybridization affinity of the sample 
of nucleic acids to nucleic acid probes that are complementary to a portion of a first allele; 

inputting a second hybridization value indicative of hybridization affinity of the 
sample of nucleic acids to nucleic acid probes that are complementary to a portion of a second 
allele; and 
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calculating an estimation of the fraction of the first and second alleles in the sample of 
nucleic acids by dividing the first hybridization value by a sum of the first and second 
hybridization values. 

5 10. The method of claim 9, further comprising subtracting from the first and 

second hybridization values a third hybridization value indicative of hybridization affinity of 
the sample of nucleic acids to nucleic acid probes that include at least one nucleotide that 
differs from both the first and second alleles. 

10 11. The method of claim 1 0, further comprising selecting the third hybridization 

value from a plurality of hybridization values, each of which is indicative of hybridization 
affinity of the sample of nucleic acids to nucleic acid probes that include at least one 
nucleotide that differs from both the first and second alleles. 

15 12. The method of claim 1 1 , wherein selecting the third hybridization value 

includes selecting a hybridization value indicative of lowest hybridization affinity to the 
sample of nucleic acids. 

1 3 . The method of claim 1 0, further comprising: 

20 calculating a first median of a plurality of hybridization values that are indicative of 

hybridization affinity of the sample of nucleic acids to nucleic acid probes that include at least 
one nucleotide that differs from the first allele; 

calculating a second median of a plurality of hybridization values that are indicative of 
hybridization affinity of the sample of nucleic acids to nucleic acid probes that include at least 
25 one nucleotide that differs from the second allele; and 

selecting the third hybridization value to be equal to the first or second median that is 
indicative of highest hybridization affinity to the sample of nucleic acids. 

14. The method of claim 1 0, wherein the subtracted first and second hybridization 
30 values are set equal to zero if they are less then zero. 
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15. The method of claim 9, further comprising: 

generating a plurality of estimations of the fraction of the first and second alleles in a 
plurality of samples of nucleic acids; and 

clustering the plurality of estimations to form a fraction pattern. 

5 

1 6. A computer program product that estimates a fraction of alleles in a sample of 
nucleic acids, comprising: 

computer code that receives a first hybridization value indicative of hybridization 
affinity of the sample of nucleic acids to nucleic acid probes that are complementary to a 
1 0 portion of a first allele; 

computer code that receives a second hybridization value indicative of hybridization 
affinity of the sample of nucleic acids to nucleic acid probes that are complementary to a 
portion of a second allele; 

computer code that calculates an estimation of the fraction of the first and second 
15 alleles in the sample of nucleic acids by dividing the first hybridization value by a sum of the 
first and second hybridization values; and 

a computer readable medium that stores the computer codes. 

17. The computer program product of claim 1 6, wherein the computer readable 
20 medium is selected from the group consisting of floppy disk, tape, flash memory, system 

memory, hard drive, and a data signal embodied in a carrier wave. 

18. A method of genotyping a sample of nucleic acids, comprising: 
generating a plurality of estimations of the fraction of first and second alleles in a 

25 plurality of samples of nucleic acids; 

clustering the plurality of estimations to form a fraction pattern; 
calculating an estimation of the fraction of the first and second alleles in a sample of 
nucleic acids; and 

genotyping the sample of nucleic acids by comparing the estimation of the fraction of 
30 the first and second alleles in the sample of nucleic acids to the fraction pattern. 
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19. The method of claim 1 8, wherein the calculating an estimation of the fraction 
includes: 

inputting a first hybridization value indicative of hybridization affinity of the sample 
of nucleic acids to nucleic acid probes that are complementary to a portion of the first allele; 
5 inputting a second hybridization value indicative of hybridization affinity of the 

sample of nucleic acids to nucleic acid probes that are complementary to a portion of the 
second allele; and 

calculating the estimation of the fraction of the first and second alleles in the sample of 
nucleic acids by dividing the first hybridization value by a sum of the first and second 
1 0 hybridization values. 

20. The method of claim 1 8, wherein genotyping the sample of nucleic acids 
includes determining if the estimation of the fraction of the first and second alleles in a sample 
of nucleic acids is within a predetermined range of a cluster in the fraction pattern. 

15 

21 . The method of claim 1 8, wherein the fraction pattern includes a cluster for a 
homozygous first allele, a heterozygous first and second alleles, and a homozygous second 
allele. 



20 22. The method of claim 1 8, further comprising identifying a single nucleotide 

polymorphism at an interrogation position based on the number of clusters in the fraction 
pattern. 

23. The method of claim 1 8, further comprising confirming that there is a single 
25 nucleotide polymorphism at an interrogation position based on the number of clusters in the 
fraction pattern. 



24. A computer program product that genotypes a sample of nucleic acids, 
comprising: 

30 computer code that generates a plurality of estimations of the fraction of first and 

second alleles in a plurality of samples of nucleic acids; 
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computer code that clusters the plurality of estimations to form a fraction pattern; 

computer code that calculates an estimation of the fraction of the first and second 
alleles in a sample of nucleic acids; 

computer code that genotypes the sample of nucleic acids by comparing the estimation 
of the fraction of the first and second alleles in the sample of nucleic acids to the fraction 
pattern; and 

a computer readable medium that stores the computer codes. 

25. The computer program product of claim 24, wherein the computer readable 
medium is selected from the group consisting of floppy disk, tape, flash memory, system 
memory, hard drive, and a data signal embodied in a carrier wave. 
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